Induced Ferromagnetism due to Superconductivity in Superconductor-Ferromagnet 

Structures. 
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We consider a superconductor- ferromagnet (S/F) structure and assume that above the super- 
conducting transition temperature Tc the magnetic moment exists only in F. In a simple model of 
the ferromagnet (the exchange field is of the ferromagnetic type for all energies) we show by an 
explicit calculation that below Tc the magnetic moment may penetrate the superconductor. In this 
model its direction in S is opposite to the magnetization of free electrons in the ferromagnet. The 
magnetization spreads over a large distance which is of the order of the superconducting coherence 
length and can much exceed the ferromagnet film thickness. At the same time the magnetic 
moment in the ferromagnet is reduced. This inverse proximity effect may explain the reduction in 
magnetization observed in recent experiments and may lead to a strong interaction between the 
ferromagnetic layers in F/S/F structures. 



Penetration of the superconducting condensate into a 
normal metal in the superconductor (S)-normal metal 
(N) heterostructures is a well established proximity ef- 
fect. The latter is a long range effect because the ampli- 
tude of the condensate decays in the normal metal very 
slowly with a characteristic length ^jv which in the dirty 
limit is equal to ^at = Djss jlixT { Dn is the diffusion 
coefRcient in the normal metal and T is the temperature). 
At low temperatures this length can be very large. At the 
same time, the order parameter A in the superconductor 
near the S/N interface is suppressed. The magnitude of 
the suppression depends on the parameters characteriz- 
ing the system such as the S/N interface transparency, 
the thickness of the S and N layers, etcQ]. The proxim- 
ity effect arises also in superconductor- ferromagnet (S/F) 
structures. While the superconducting condensate con- 
sists of paired electrons with opposite spins, the exchange 
field J in the ferromagnet tends to align them and break 
the Cooper pairs. The penetration length of the con- 
densate into the ferromagnet is usually much smaller 
than and in the dirty limit is equal to = \/DfJJ 
{Dp is the diffusion coefRcient in the ferromagnet). Since 
the exchange energy J is much larger than T, we come to 
the inequality <^ (h^ the clean limit when r J >> 1, 
the penetration length ^p is determined by the mean free 
path / ~ vpT, where vp is the Fermi velocity[illl3). A 
strong exchange field suppresses also the superconduct- 
ing order parameter A in the superconductor. 

The situation changes when the magnetization M in 
the ferromagnet is not homogeneous. In this case a triplet 
component of the condensate with a non-zero spin pro- 
jection arises and penetrates into the F- region over a 



long distance of the order Dp 12-kT (see Refs. 0,0|). 
The effect of the penetration of the superconducting con- 
densate into the ferromagnet and the suppression of the 
superconductivity (a decrease of the critical temperature 
Tc of the superconducting transition) in S/F structures, 
i.e. the proximity effect, has been the subject of many 
works during the last decades (see for example (St iSfc ISl 
and the review IQj for more references). 

So, the penetration of the superconductivity into the 
normal metal or ferromagnet is by now a very well studied 
phenomenon. However, one can ask the same question 
about the ferromagnetism: Can the ferromagnetic order 
penetrate the normal metal or superconductor over long 
distances? Surprisingly, this question has hardly been 
addressed. Some indications of the effect can be found 
in numerical works Refs. 0,^3- In those works only the 
density of states for each spin direction as a function of 
the energy was presented, however the magnetization was 
not calculated. In addition the induced magnetization 
( "magnetization leakage") was calculated in Ref. 
However the results obtained in the latter paper gener- 
ally speaking differ drastically from ours. They found a 
" magnetization leakage" , that is the magnetic moment 
of free electrons Me spreads into the S region over a dis- 
tance of the order ^5 changing its sign at some distance 
from the S/F interface. 

We consider a simple model assuming a mean field 
approximation for the ferromagnet and superconductor. 
The mean field order parameter in S is the energy gap 
A, and in F it is the exchange field J which is assumed 
to be of the ferromagnetic type and small compared with 
the Fermi energy. In different limiting cases where ana- 
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lytical formulae can be obtained we find completely dif- 
ferent behavior: For temperatures below Tc, the mag- 
netization of free electrons in the F layer decreases 
and the induced magnetization in the S region is negative 
(that is, the magnetization variation has the same nega- 
tive sign in both regions). Our analytical considerations 
show that no change of sign of the induced magnetization 
takes place. This behavior is in agreement with the re- 
duction of magnetization observed in the experiments of 
Refs.jl^[l3| and can be explained by the simple physical 
picture we present below. 

In the case of F/N systems the ferromagnetic order- 
ing penetrates over short distances since the exchange 
interaction is local. In this paper we show that the sit- 
uation may be different for S/ F structures and present 
arguments that the magnetic moment can penetrate into 
the superconductor over long distances of the order of 
the superconducting coherence length. This effect can 
be called the inverse proximity effect. The reason why 
the magnetic moment aligned in the direction opposite to 
the magnetization in the F film Mg penetrates the super- 
conductor can rather easily be understood qualitatively. 
This effect is due to the fact that the Cooper pairs have 
a large size of the order of — \J Dsl'^i^Tc- Suppose 
that the F layer is thin (see inset of Fig. ^ and let us 
assume that the Cooper pairs are rigid objects consist- 
ing of electrons with opposite spins, such that the total 
magnetic moment of a pair is equal to zero. Of course, 
the exchange field should not be very strong, otherwise 
the pairs would break down. It is clear from this simple 
picture that pairs located entirely in the superconduc- 
tor cannot contribute to the magnetic moment of the 
superconductor because their magnetic moment is sim- 
ply zero. Nevertheless, some pairs are located in space 
in a more complicated manner: one of the electrons of 
the pair is in the superconductor, while the other moves 
in the ferromagnet. These are those pairs that create 
the magnetic moment in the superconductor. The direc- 
tion along the magnetic moment Mg in the ferromagnet 
is preferable for the electron located in the ferromagnet 
and this makes the spin of the other electron of the pair 
be antiparallel to Mg. This means that all such pairs 
equally contribute to the magnetic moment in the bulk 
of the superconductor. As a result, a ferromagnetic or- 
der is created in the superconductor and the direction of 
the magnetic moment in this region is opposite to the 
direction of the magnetic moment Me in the ferromag- 
net. Moreover, the induced magnetic moment penetrates 
over the size of the Cooper pairs ^s- From this point of 
view it is difhcult to understand the numerical results of 
Ref. where the induced magnetization in the S re- 
gion near the S/F interface has the same sign as in the 
ferromagnet. The magnetic moment in the ferromagnet 
is decreased because the density-of-states in F is reduced 
due to the proximity effect. This occurs in a way similar 
to a suppression of the Pauli paramagnetism in super- 




FIG. 1: Spatial dependence of the magnetization in the whole 
sy stem . Here ^whs = 0.5, 7f = 7f/Co = 0.1 (^o = 
^ Ds I'^Tc), J/Tc = 15 and dF/Co ~ 1- Inset: Schematic view 
of the inverse proximity effect in a S/F system (for discussion 
see text). 

conductors (the exchange field plays the role of a strong 
magnetic field acting on spins). At the same time the 
concentrations of free electrons with spin up and down 
in F remain unchanged; when we are saying about the 
penetration of Cooper pairs into the ferromagnet F, we 
mean that superconducting correlations are established 
in F due to the proximity effect. 

Having presented the qualitative picture, we calculate 
now the magnetization variation due to free electrons (the 
conduction band electrons) 

6M, = ^iB5NM (1) 

below Tc in both layers of the S/F system shown in the 
inset of Fig^ Here fXB is an effective Bohr magneton. 
We assume that the magnetic moment Mg is parallel to 
the interface as it takes place in the experiment 0, Q| 
and is homogeneous in the F layer. As we have found 
previously, Ref.^^, in this case only the singlet compo- 
nent and the triplet one with the zero spin projection 
on the direction of Me exist in the system. Both com- 
ponents penetrate into the ferromagnet over the short 
distance • If the S/F interface transparency is low or 
the conductivity of the S film is much higher than the 
conductivity of the F film, the suppression of the order 
parameter A is not essential and the superconducting 
properties remain almost unchanged. 

The quantity 5Nm can be expressed in terms of the 
quasiclassical normal Green function g 

p UJ— — 00 

(2) 

where v = ppm/ is the density of states at the 
Fermi level, (73 is the third Pauli matrix and uj — 
7TT{2n + 1) is the Matsubara frequency. The normal 
Green function g is a matrix in the spin space. In the 
considered case of an uniform magnetization it has the 
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form g = go ■ <Jq + gz ■ cr^- This matrix is related to 
the Gor'kov anomalous matrix Green function / via the 
normalization condition 



1 



(3) 



The matrix / describes the superconducting condensate. 
In order to visualize how our results are obtained, we 
consider first the simplest case when the condensate func- 
tion / is small in F and is close to its bulk value in the 
superconductor. We analyze the dirty case when the 
Usadel equation can be applied. This means that s-wave 
superconductors are considered. These are described by 
the standard BCS Hamiltonian with account for the ex- 
change field (in the ferromagnet) acting on the spins of 
the free electrons. We write the Hamiltonian in the form 



H — Ho + Hs + Hp 



(4) 



Here Ho is the one-particle Hamiltonian which includes 
the impurity scattering term Uimp 

Ho='^ {a% [ip^pp' + Uimp6ss']as'p'] (5) 

where = /2m — ej- is the kinetic energy counted 
from the Fermi energy ep- The second term in Eq.Q) is 
the standard BCS Hamiltonian for the superconductor 
written in the mean-field approximation 



Hs 



E {^4p^tp + c.c] 



(6) 



where s = — s (s = ±1) and p = —p. The last term in 
Eq.Q) describes the ferromagnetic interaction in F. We 
accept the simplest form of this part of the Hamiltonian 
(the mean-field approximation) 



where S = J2a ^a^i^ ~ '^a. is the spin of a partic- 
ular ion. A constant Ji is related to J via the equation: 
J = Jin M So , where njvfis the concentration of mag- 
netic ions and Sq is a maximum value of S n (w e consider 
these spins as classical vectors; see Ref. lig). In this 
case the magnetization is a sum: M = M/oc+Mg, and 
the magnetization Mg may be aligned parallel ( Ji > 0, 
the ferromagnetic type of the exchange field) to M or 
antiparallel ( Ji < 0, the antiferromagnetic type of the ex- 
change field). In the following we will assume a ferromag- 
netic exchange interaction (Me and M are oriented in the 
same direction). The case of antiferromagnetic coupling 
will be briefly discussed below. In principle one can add 
to Eq. (jHJ the term ^| {Sa * Sb} which describes a 
direct interaction between localized magnetic moments 
[itIIi^ . but this term does not affect final results. 

Starting from the Hamiltonian and using a stan- 
dard approach 18], one can derive the Usadel equation. 
In the case of a low S /F interface transparency this equa- 
tion can be linearized. Then, the function / is obtained 
from the hnearized Usadel equation (see e.g. Ref. 15]) 



0, in the F layer 



and 



dlx^fs - «s'5/s = ^(2^)0-3, in the S layer. 



(9) 



(10) 



Here k\ = 2{\uj\ ^ iJsgnu)/ Dp, K5 = 2^/L0^ + ^?/Ds 
and 5fs is a deviation of the function fs from its bulk 
(BCS) value Jbcs, i-c 5fs = fs - fscs, fncs = 
Ibcs ■ o's. The functions fp± are the elements (1,1) and 
(2,2) of the matrix fp. The function JC{x) contains the 
correction 5/S.{x) to the order parameter A. This term 
is not relevant in our calculations since only the com- 
ponent of / proportional to (Tq contributes to the mag- 
netization (see below, Ea. ltT7|l 'l. Eqs. H9I1()|I should be 
complemented by the boundary conditions that can be 
written for small fs.p as 



Hp = - 



{p-A 



J {a 



sp 



n * (Ts 



(7) 



where the exchange energy is assumed to be positive 
for all energies (the ferromagnetic type of interaction), 
n is the unit vector parallel to the magnetization of the 
ferromagnet. The magnetization M of the ferromagnet 
is proportional to the exchange energy J. If the contri- 
bution of free electrons strongly dominates (an itinerant 
ferromagnet), one has M = Me. 

If the polarization of the conduction electrons is due 
to the interaction with localized magnetic moments, the 
Hamiltonian Hp may be written in the form |16L llTl Il9| 



Hp 



-Ji {a+ S : 



'P' } 



(8) 



dxSfs = {9 BCS ■ If - gBCsfBCsSgnuj • o-3)/75 (11) 
dxfp = -{^/lF)fs , (12) 

where js,f — Rb^s.F, Rb is the S/F interface resistance 
per unit area, as^p is the conductivity of the S or F re- 
gion, and gBCS = gBCS • ^o- The BCS functions have 
the well known form (see for example [l^ ) 



gBCS = Lu/Vuj^ + A^, Ibcs = A/iy/iu^ + (13) 

The matrix function of the superconducting conden- 
sate / can be represented in the form 



/ = /3<5'3 + /oo'o 



(14) 



for both regions. The component describes the singlet 
condensate, whereas fo stands for the triplet component 



4 



with the zero projection of the total spin of the pair on 
the direction of the magnetic moment M. This functions 
are related to /± through: fo,3{x) = (l/2)(/+(a:)±/_ (.t)) 
The other components of the triplet condensate arise only 
if M in the ferromagnet is inhomogeneous. 

Solving Eas. (|9ll0|l with the boundary conditions, Eqs. 
(tTT|l and l(T^ . we find easily 

fF±{x) = b±exp(-K±x) (15) 
fso{x) = -aoexpinsx) . (16) 

Here b± = ±fBCs/ilFK±) and oq = 
9BCsfFo{0)/{'ysKs)- As follows from Eqs. ^ 
and H16(l the functions /i?± and fso{x) are small pro- 
vided that Rp/Rb « 1 and {Rp / Rb){Rs / Rb) « 1, 
where Rf,s — ^f,s/o'f,s are the resistances (per unit 
area) of the F(S) of lengths ^f,s- In order to calcu- 
late the magnetization we have to find the function 
33 ~ Tr&3g/2 (see Eq. (jSJ). The latter is related to 
the functions /0.3 in the F and S region through the 
normalization condition, Eq.|(31l, and is given by 

gF3 = fFofFsSgnw, gs3 = fBCsSfso/gBSC ■ (17) 

As it has been discussed in Ref. 0] , the functions /so and 
fpo corresponding to the triplet component of the con- 
densate are odd function of uj while the singlet compo- 
nents fscs and fp3 are even functions. Thus, according 
to Eas. ljlTII the functions gF3 and gs3 are even functions 
of u) ( gBSC is odd in w). This means that the sum over 
the frequencies in Eq.Q is not zero and the proximity 
effect leads to a change 6Me of the magnetization in both 
F and S layers (above Tc the magnetization in S is zero). 

After the qualitative discussion we have come to the 
conclusion that the net magnetization due to the inverse 
proximity effect must be negative. The explicit calcula- 
tion based on Eqs.© and H15|l - (|17|l confirms this result 
which is shown in Fig. ^for some values of the parame- 
ters. We see that SMg is negative, i.e the magnetization 
of the ferromagnet is reduced and the superconductor 
acquires a finite magnetization in the opposite direction. 
The change of the magnetization 5M{x) extends over the 
length Kj^, which may be much larger than the thickness 
of the F layer. This effect is another manifestation of the 
existence of the triplet component of /. 

The inverse proximity effect considered here may be 
relevant for several experiments on measurements of the 
magnetization in the S/ F structures 0, Q| . In these 
experiments it was found that the magnetization started 
to decrease when crossing the superconducting critical 
temperature Tc from above. The authors of these exper- 
iments compared the data with the theoretical results of 
Refs. |l9t |2(]j that were obtained under the assumption 
that the ferromagnetic order in the F thin layer might 
be modified due to the proximity effect leading to the 
so called cryptoferromagnetic state. In these works only 



the contribution of localized moments to the magnetiza- 
tion was taken into account. Our calculations show that 
the conduction electrons can give an additional contribu- 
tion. The inverse proximity effect leads to an additional 
reduction of the magnetization M^, and may serve as an 
alternative explanation for the reduction of the magneti- 
zation observed experimentally 0, 0| . 

Let us analyze now an interesting case that may be 
relevant to the experimental situation of Ref. 1^. We 
assume that the thickness of the F layer dp is small com- 
pared to ^p and that the Green functions gs and /s are 
close to the bulk values gBSC and Jbcs- The latter as- 
sumptions is valid if the coefficient ^f/^s — <^fI<^s is 
small enough. In this case all functions in the F region 
are not necessarily small but they are almost constant 
in space. Therefore we can average the exact Usadel 
equation over x taking into account exact boundary con- 
ditions. Proceeding in this way, we get for the diagonal 
elements g± and /± of the matrices g and / 

gp± = i^±liuj±, !f± = ±^bFfBCs/Cuj± (18) 



where u)±=oj + ei^pgBCsTiJ , C± ^ \J^±- {^bFfBCsY, 
CbF = Df / i'^jFdp)- One can see that in the limiting 
cases of small and large energy ebF the functions gF±, 
fpj. describe a superconducting state with the energy 
gap equal to ebp if e&_F << A ( a subgap in the excita- 
tion spectrum) and to A in the opposite case. In both 
cases the position of the energy gap is shifted with re- 
spect to e = (the Matsubara frequencies are related to 
e via uj — —ie). It can be easily shown that the func- 
tion gp3 that determines the magnetization, Eas. (|2ll7|l . 
equals zero for ebP = (very small S/F interface trans- 
parency) and for very large values of ebF (a perfect S/F 
contact). This dependence of SMp on Rb leads to a non- 
monotonic behavior of the change of the magnetization 
SMp,s- In FigI21we show the temperature dependence of 
SM p cjO d) for values of the parameters similar to those of 
Ref. [iJI . We see that the decrease of the magnetization 
may be of the order of 10% and larger. This results corre- 
late with the experimental data of Ref. 0| (see footnote 
[^). We have checked that |/so| < 1 for the parameters 
in Figia 

We also present here analytical formulae for the ra- 
tio rs,F = dMs,p{0)/Mpo using Eq. ((TH|l and consid- 
ering the case of low temperatures {T << A); SMs^f 
are the magnetization variations in the S and F films 
and MpQ is the magnetization in the F film above Tc- 
The relation between MpQ and J depends on a particular 
model of the ferromagnet. For example, in the simplest 
model of the ferromagnet with a constant and positive J 
we have for an itinerant ferromagnet MpQ= g^B^'J (see 
[2ll|V To simplify the expressions for rs,F we assume also 
that J « ebF ~ {Dp/(fp){Rp/Rb)idF/^F) (this limit 
may correspond to the experiment [14] l. In this case we 
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8M/M, 




-0.2 , 



T/Tc 



FIG. 2: Temperature dependence of SMs/Mfq (lower curves) 
and 5Mf /Mfo (upper curves) for the following values of 7f = 
Jf/^o (eo = y/Ds/2Tc): 7f = 0.1 (solid line), -yj. = 0.3 
(dashed line), and jf ~ 0.5 (dot-dashed line). Here 7f/7s ~ 
0.5, J/Tc = 20, dF/^o = 0.1 and A = 3 (see [18]) 



obtain 



rs « -1.67^A4/Ds 
rp w — 7rA/2ef,i? 



(19) 
(20) 



For estimations of the parameters one can take experi- 
mental values from Ref . where a " weak" ferromagnet 
Pd(^i-x)F^x was used. One gets Dp /dp — lOOOX for 
dp = 20A. The Curie temperature which may be of the 
order of J varied from 90 to 250 K. The barrier (interface) 
resistance Ri, is not known, but one can give a crude esti- 
mation noting that (Rp / Rb)idp /^p) w Ttj., where T^^ is 
the transmission coefficient which varies from very small 
values to a value of the order 1. 

In conclusion, we have demonstrated the existence of 
the inverse proximity effect in S/F structures. Due to 
the presence of the superconductor the magnetization in 
the ferromagnet with the ferromagnetic type of the ex- 
change interaction is reduced and a magnetic moment is 
induced in the superconductor below Tc- Its direction 
is opposite to the direction of the magnetic moment in 
the ferromagnet and spreads over the superconducting 
coherence length ^g. This distance can be much larger 
than the F film thickness. The effect discussed may be 
the reason for a reduced magnetization observed in S/F 
structures leading to a frequency shift of the magnetic 
resonance ^1 4] . This conclusion is changed in the case of 
ferromagnets with the antiferromagnetic interaction be- 
tween free electrons and localized moments (negative Ji 
in Eq.JSl). If a contribution of localized moments to 
the total magnetization in such ferromagnets dominates 
(M/oc > Mf.), the magnetisation Mg is opposite to M 



and therefore the induced magnetization variation in the 
superconductor JMg will be parallel to M. 

Note: After completing this work we became aware of 
the paper Ref.|2^ where the magnetization leakage into 
the S layer was numerically calculated for a ballistic S/F 
structure. In this case the magnetization penetrates the S 
layer over distances of the order of the Fermi wave length. 
We are not interested in small scales of this order. 

We would like to thank SFB 491 for financial support. 
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